On the occurrence of interference in special two-slit experiments 
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We investigate numerically the occurrence of the interference fringes in experiments where an 
initial Gaussian wave packet evolves inside billiard domains of several shapes with two narrow slits on 
a side. Contrarily to what claimed in the literature, our results seem to show that the occurrence of 
interference fringes does not depend on the billiard integrability but on a spatial reflection symmetry 
concerning both the billiard domain and the initial Gaussian wave packet. Indeed, whether the 
billiard is regular or chaotic, we find clear interference fringes when this symmetry is verified whereas 
when it is not verified, we find that the interference patterns are perturbed in various measure 
according to the effects of the symmetry violation. 

PACS numbers: 03.65.Ta, 03.65.Yz, 05.45.Mt 



I. INTRODUCTION 

In a recent work lj, Casati and Prosen have shown 
an interesting numerical simulation concerning the 
quantum-mechanical time-evolution of an initial Gaus- 
sian wave packet inside a billiard domain having the fol- 
lowing shapes: 

i) isosceles right triangle 

ii) as above but with the hypotenuse replaced by a 
circular arc. 

In both cases the domain has two narrow slits on a 
cathetus, through which the initial wave packet leaks out 
bit by bit. The set-up of the experiment is intriguing, 
because it involves properties of classical dynamics and 
properties of quantum dynamics. The former are repre- 
sented by the shapes of the billiards, where the classi- 
cal ray dynamics can be regular (case i) or fully chaotic 
(case ii)), whereas the latter are represented by the in- 
terference of the two wave packets coming out from the 
slits. The experiment gives thus the opportunity to inves- 
tigate two questions which, although debated for a long 
time, still remain far from being fully understood, i.e. the 
manifestations of classical chaos in quantum mechanics 
(quantum chaos) and whether an external decoherence 
mechanism is necessary in order that a system presents 
the quantum to classical transition 0, Q . The result in 
is that interferences fringes occur only in the case of 
the integrable billiard i). Casati and Prosen have drawn 
thus the conclusion that their experiment shows: 



1) "a vivid and fundamental illustration of the mani- 
festation of classical chaos in quantum mechanics" 

and that 

2) "considering a pure quantum state in absence of 
any external decoherence mechanism, internal dy- 
namical chaos can provide the required randomiza- 
tion to ensure quantum to classical transition in the 
scmiclassical region" . 

An investigation similar to 1] has been carried out p| 
in an Aharonov-Bohm two-slit set-up by a semiclassical 
calculation of the outgoing probability density current. 
This further investigation substantially confirms the fact 
that the occurrence of interference fringes depends on 
the billiard integrability but it is critical of conclusion 
2). Impressed by the conclusions 1) and 2), in the present 
paper we report on a numerical simulation similar to that 
in [lj but concerning a much larger family of billiards. 



II. DESCRIPTION OF THE NUMERICAL 
EXPERIMENT 

A billiard B (in Fig. 1 represented by a square) is set 
inside a rectangular integration region R. The borders 
of B arc defined by a thin potential barrier Vs(x,y), 
sufficiently high to reduce as much as possible the tunnel 
effect. Two slits a and b on the lower side I of B are 
obtained setting Vb(x, y) = along two small intervals 
of I. The numerical simulation consists in the integration 
of the Cauchy initial value problem for the Schrodinger 
equation [jj 
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i^ix.V.t) = H^{x,y;t) 
ip(x,y,0) = ipo{x,y), 
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and tpo is the normalized Gaussian wave packet 



The intensity I{x) on the screen S at x, i.e. the prob- 
ability density gone across the point x of S during the 
experiment, is then evaluated by 



I ( x ) = / j y (x,y;t)dt, 
Jo 



(4) 



where T is the duration of the experiment, j y the com- 
ponent along the y-axis of the probability density cur- 
rent and y the position of the screen. The integration of 
the problem (1) has been carried out through the eval- 
uation of the time-evolution operator e~ lHt by a finite 
difference method, where e~ lHt has been approximated 
by the exponential series truncated at the fourth order. 
The order of this truncation is crucial for the accuracy of 
the method. Indeed, the truncation at the first order em- 
ployed in previous papers gives rise to serious problems 
of stability and probability conservation (see Ref. and 
references therein). Strictly speaking the integration re- 
gion of (1) is the whole space M 2 , but in our calculations 
we have reduced this space to the finite rectangular region 
R. This reduction causes reflections from the border of 
R. In order to avoid this problem, we have employed the 
same methods in that is damping the wave function 
by means of an imaginary potential — iVa{x, y), where 
VA(x,y) > inside an "absorbing layer" A (see Fig. 1) 
along the border of R and zero elsewhere. The integra- 
tion of (1) with this method turns out to be very accurate: 
the total probability is conserved and the transmission 
through the boundary of B and the reflection from the 
sides of R are negligible. It should be observed that a way 
alternative to (1) to descrive the time development of tpo 
is that of considering a Cauchy initial value problem like 
(1), but relative to the free-Hamiltonian — ^A, with in 
addition the boundary condition ip(x,y-t) — 0, Vt > 0, 
on the border of -B^(|- The method (1), employed in the 
present paper, has a twofold advantage. The first one 
is of a formal nature. Whatever the billiard shape may 
be, the description of the time-evolution of "00 is given 
by the problem (1), which is defined in the whole space 
R 2 and presents a Hamiltonian operator which is always 
(see Ref. 0, p. 182) self-adjoint in the space L 2 (R 2 ). 
On the contrary, in the other method one should prove 
the self-adjointness of — in each billiard domain with 
the boundary condition tp — on its border. The second 
advantage of the method (1) is of a practical nature: one 
can handle any billiard shape without any difficulty In- 
deed, there is no need to impose each time the tiresome 
boundary condition ip = on the border of the billiard 
and all that changes with the billiard shape is only the 
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FIG. 1: Experiment set-up. H height and L length of the 
rectangular integration region, A absorbing layer, B billiard, 
I lower side of B, O origin of the coordinate system, a and b 
slits, S screen, G gaussian wave packet (3), k = (k x , k y ) wave 
vector of G. 



Parameter 


Value 


H 


1.6 


L 


1.2 


length of / 


1 


width of the potential barrier Vb(x,u) 


0,008 


height of the potential barrier Vb(x,j/) 


10 b 


w (slit width) 


0.012 


d (slit distance) 


0.1 


s (screen distance from I) 


0.3 


S (space step size) 


0.002 


t (time step size) 


10~ b 


ro (center of the packet (3)) 


(0.0,-0.25) 


(To (rms width of the packet (3) ) 


0.09 


| fe|| (magnitude of k ) 


180 



TABLE I: Parameters employed in our numerical simulation. 



expression for VB(x,y) in (1). 

As regards the parameters, we have chosen values de- 
liberately close or equal to those in [|. See Table I and 
also Fig. 1. 

The choice of the values for ctq and ||fc|| yields a 
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FIG. 2: Specimens of our results. Column I: Billiard with the initial wave packet (3) inside. The center of (3) has coordinates 
xo = 0, yo = —0.25. The arrow denotes the wave vector k in (3). The components of k are k x = 0, k y — 180 in a),b),e), 
k x = —113, k y — 140 in c), d) and k x = 113, k y — 140 in /). The center of the ring in the billiards has coordinates x = 0, y — —0.6 
in b),d) and x = 0.2, y — —0.7 in e),/). Column II: Intensity (4) as a function of the position on the screen. Column III: 
Sample plot of cos(A<^(t)) as function of time, where Aip(t) is the phase difference of the solution ip(x,y;t) of (1) at the slits. 
In each row the figures in Columns II, III are relative to the billiard in Column I. 
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FIG. 3: Wave fronts of the solution ip(x,y;t) of problem (1) 
in the case /). Left half: one slit open, right half: both slits 
open. Quite similar figures are obtained for the other billiards. 



Gaussian wave packet of good monochromaticity, indeed 
~ 3.1 • 10~ 2 , where Ok is the rms width of (3) in 
the momentum space . The values of d and w have been 
chosen consequently in order to ensure a good visibility 
of the fringes. 



III. RESULTS AND DISCUSSION 

In Fig. 2 we report sample results of our numerical 
experiment, which has been stopped when the probability 
density leaked out from the slits was of order 85 per cent 
of the initial probability inside B ^lj ■ 

The results a) and b) show that the interference fringes 
occur whether when the billiard is regular (square bil- 
liard a)) or when it is chaotic (Sinai billiard 6)), provided 
that the wave vector fc of the initial wave packet (3) is 
directed along the y-axis^3- When fc is varied of the 
same finite quantity, both in the square billiard c) and 
in the Sinai billiard d), the situation, independently of 
the value of this variation, changes drastically: the inter- 
ference fringes almost disappear and are replaced by two 
comparable intensity patterns. The plots in Column III 
point out clearly the new state of things. Thus, contrar- 
ily to what claimed in 1] and confirmed in we can 
say that the occurrence of the interference fringes does 
not seem to be correlated with the integrability of the 
billiard but with the direction of the wave vector fc. On 
the other hand, the cases e),/) show that the only di- 
rection of fc is not sufficient to determine the appearance 
of the fringes, indeed, in these latter cases, whatever the 
direction of fc may be, the interference fringes never oc- 
cur and quite similar intensity patterns are obtained on 
varying the direction of fc 01 • Taking all these facts into 
account, we can now find out what actually determines 
the occurrence of the fringes. Indeed, one expects that 
these appear when the time development of the initial 



Gaussian wave packet (3) gives rise at the slits a, b to 
two respective coherent waves. Seen the symmetric posi- 
tion of the slits with respect to the y-axis (see Fig. 1), a 
sufficient condition for this coherence is that the solution 
ip(x,y;t) of the Cauchy problem (1) is such that 

i>(x, D\ t) = ip(-x, y; t), Vi > 0, Mx, y e R. (5) 

We now see that this property is verified if the potential 
barrier Vb and the initial Gaussian wave packet (3) are 
invariant for a spatial reflection w.r.t. the y-axis, i. e. 

Vb(x,v) = V B (-x,y), a,ndip (x,y) = 4> (-x,y) Vx,y e B. 

(6) 

Indeed, performing such a transformation, the Cauchy 
problem (1) becomes 

f i%il>(-x,v,t) = [-±A + V B {-x,y;t)}i)(-x,y;t) 
1 ip(-x,V,0) = ip (-x,y), 

(7) 

and hence, if condition (6) is satisfied, we have (5). 
Henceforth, we shall call (6) symmetry condition (SC). 
The results a) -f- /) can be now explained as follows. 

We regard the slits a, b as point sources of two re- 
spective wave packets ij) a (x,y;t), , iljb(x,y;t) such that 
i/j(x,y;t) = ip a (x,y;t) + ip b (x,y;t), V (x,y) belonging 
to the region below the slits, V t > 0, and ip(x, y; t) = 
if> a (x,y;t),ip(x,y;t) = ip b (x,y;t), Vi > 0, at the slits. 
When the SC is verified (cases a), b)), the relation (5), 
as we said, holds. Thus, in particular, ip a (x,y;t) = 
tpb{x,y;t), Vi > 0, at the slits. In other words, the slits 
radiate two coherent waves (see the plots of the phase 
difference relative to a),b) in Column III of Fig. 2). If 
the SC is not verified, as in the cases c),d),e), f), even 
if only one of the two conditions in (6) is not satisfied, 
the relation (5) does not hold any more and we do not 
have any guarantee that ip a (x,y;t) = ip b (x,y,t) at the 
slits and hence that the slits radiate two coherent waves. 
Actually (see the plots of the phase difference relative to 
c), d), e), /) in Column III of Fig. 2), the two waves emit- 
ted do not show now any definite phase relationship with 
each other. In any case (SC satisfied or not satisfied), 
the waves ">p a ,ipb radiated by the slits are circular (see 
Fig. 3) and have a wave length A which can be estimated 
numerically of order 3 • 10~ 2 . It is now fairly easy to ex- 
plain the appearance of the interference fringes. Let P 
be a point on the screen S and let r a , r b be the respective 
distances from the center (—d/2, 0) of the slit a and from 
the center (d/2,0) of the slit b (see Fig. 1). Then, the 
phase difference Aip(P, t) of ip a and ipb at P and at time 
t can be written as 

A^(P; t) = <p a (-d/2, 0; t) - <p b (d/2, 0; t) + • 2vr, 

(8) 

where <p a (— d/2, 0; t) and ip b (d/2,Q;t) are, respectively, 
the phases of ip a and tpb at the slits. Now, as we have said, 
when the SC is verified <p a (-d/2, 0; t) - fb(d/2, 0; t) = 
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at any t, and thus by (8) the interference fringes occur, 
independently of the the billiard shape. 

When the SC is not verified, the difference 
ip a (—d/2, 0; t) — (fb(d/2, 0; t) is a random function of time 
and thus, whatever the billiard shape may be, the fringes 
are perturbed to a different extent, as we shall discuss 
better below. Contrarily to the conclusion 1) in paper 
, we thus are led to say that the occurrence of the in- 
terference fringes is not determined by the integrability 
of the billiard but by the fulfillment of the SC (6) . If this 
is the case, we note that simple symmetry considerations 
permit to derive some general properties of the interfer- 
ence patterns. Indeed, we have ip a (x, V) t) = ipb{—x, y; t), 
at any t > in the region below I. Therefore at any point 
P = (x, y) on the screen 

A(p(x, y; t) = ip a (x, y; t) - <p b (x, y; t) 

ip b (-x, y; t) - (fia(-x, y; t) = -A<p(-x, y;t), V t > 0, 

and thus the interference pattern is symmetric w.r.t. the 
y-axis and the center of the screen is a maximum. In 
Fig. 4 we show typical results for the regular triangular 
billiard considered in The set-up g), except for the 
direction of fc and for the center of the wave packet, is 
just that in This set-up violates the SC whatever the 
direction of fc may be, because Vs(x,y) ^ Vb(— %, y)- 
Thus, according to what we said above, this set-up 
cannot ever give rise to a clear interference pattern. In 
fact, by varying the direction of fc we have always found 
interference patterns perturbed in various measure. 
To this family ofpatters belongs, indeed, the same 
pattern shown in [lj, because it presents minima which 
are visibly different from zero. A perfect interference 
pattern, in the sense that the minima are practically 



zero, is instead obtained in the set-up h), the only 
one (together with the other where fc has the opposite 
direction) that satisfies the SC. Notice how the plots of 
the phase difference at the slits (Column III in Fig. 4) 
show clearly that the coherence at the slits takes place 
only in the set-up h). 



As we have seen, when the SC is not verified, the inten- 
sity patterns so far shown resemble variously perturbed 
interference patterns up to almost to coincide, as in the 
cases c),d), /), with the pattern relative to two incoher- 
ent waves ip a (x, y; t), ipb(%, y, t). This variety of intensity 
patterns, according to our interpretation of the results, 
is due to the different way how the SC breaking acts on 
the decoherence of the wave packets ip a ,ipb at the slits. 
It is not easy to find out a precise relation between the 
extent of the SC breaking and the corresponding pertur- 
bation of the interference patterns. We confine ourselves 
to discuss here only the effects due to some weak SC 
breakings. Recalling (6) and observing that the condi- 
tion tp(,(x,y) = ipa(—x,y), Vx, y G B, there appearing is 
equivalent to k x = 0, xq — 0, we can say that the SC is 
violated if at least one of the following conditions holds: 
k x ^ 0,x ^ in ?/>o and V B {x,y) ^ V B (-x,y). For the 
sake of brevity, we suppose (as we have always done in 
our calculation) xq = 0, and treat in detail the effects 
due separately to small variations of fc and of the po- 
tential barrier. In the first case, let ipo and ipo be two 
initial wave packets, characterized, respectively, by wave 
vectors fc and fc. If one employs a first-order expansion of 
the exponentials exp[ik ■ r] and exp[ik ■ r] appearing, re- 
spectively, in ipo and ipo, the effects on the time-evolution 
due to the replacement of ipo with ipo can be estimated 
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FIG. 5: Results about the influence of the variations of tf>o and of the potential barrier on the interference pattern. Columns 
and rows the same as in Fig. 2. The wave vector k has components k x — —2, k y = 179.99 in i) and k x = 0, k y — 180 in l),m). 
The coordinates of the center of the ring in the billiard I) are x = 0.01, y — —0.6, whereas those of the ring in the billiard 
m) are x = 0.2, j/ = —0.7. The interference pattern in m) is relative to the probability piled up in the integration interval 
[0, 1.3 • 10 4 t]. The emission of probability in the interval [2.8 • 10 3 r, 1.3 ■ 10 4 r] is negligible. 



by 

\\e- im ^ -e- iH %\\ 2 = l^o-^oll 2 ^ 

— (kx kx) ^4 + (^3/ ^3/) B ~l~ k%) {ky ky)G, (9) 

where < A < 0.25,0 < B < 1,0 < \C\ < 0.5 and 
|| • || denotes the L 2 -norm. In virtue of (9), we can say 
that small variations of k in the initial wave packet yield 
weak effects on the interference patterns. An example of 
such an effect is shown by the case i) in Fig. 5. This 
case is obtained from a) replacing k x = 0, k y = 180 with 
k x = —2,k y = 179.99. The effect on the interference pat- 
tern is indeed barely visible and the coherence between 
ipai^Pb at the slits is, on average, partially maintained. 
We now consider the consequences on the time-evolution 
of an initial wave packet tpo due to a replacement of a po- 
tential barrier, say Vb(x, y), with another, say Vb- These 
consequences can be estimated (see Ref. @, p. 495) at 
the first order in (Vb — Vb) by 

e- ii? Vo-e- iff Vo ~ - f e- iH ^- s \VB-V B )e- iH ^^ds. 
Jo 



Contrarily to the previous case, the L 2 -norm of this quan- 
tity at each given t is not, in general, infinitesimal for 
infinitesimal variations of the potential barrier, because 
there are infinitesimal variations of the potential bar- 
rier, which do not give rise to an infinitesimal operator 
Vb — Vb- This fact is illustrated in Fig. 5 by the case I), 
which is obtained displacing the center x = 0, y = —0.6 
of the ring in the Sinai billiard d) into the new center 
x = 0.01, y = —0.6. Although this displacement is very 
tiny, the effects on the interference pattern are now quite 
evident. 

Finally, we show the case m) where, although the SC 
is broken to a high degree 0, the intensity plot, limited 
to the integration time [0, 1.3-10 4 r], reproduces a perfect 
interference pattern. This fact is not a peculiarity of 
the billiard m), but it happens similarly for any billiard, 
whatever the degree of the SC violation may be. Its 
analysis illustrates in a way particularly instructive 
that the causes of the perturbation of the interference 
fringes are the previous discussed asymmetries of the 
initial conditions and of the Hamiltonian in the problem 
(1). Indeed, during the first impact of the wave packet 
ijj(x,y;t) solution of (1) with the lower side I of the 
billiard (impact which lasts roughly 2.8 • 10 3 r), the slits 
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FIG. 6: Comparison of the intensity plots (dotted curves) at the screen in the experiments c) and /) with the plots (solid 
curves) sum of the intensities produced in two one-slit experiments, derived from the original two-slit set-up by closing the slits 
a and b alternatively. 



radiate two coherent wave packets tl> a ,ipb (see the plot 
of the phase difference in Column III of Fig. 5), because 
the time-evolution of ip(x,y;t) is not yet under the 
influence of the effects due to the direction of k and to 
the asymmetry in the billiard shape. Consequently, the 
wave fronts of ip(x,y;t) are basically the same as those 
in the initial Gaussian wave packet (3), i.e. straight lines 
parallel to I. In the time interval [2.8-10 3 r, 1.1 -lOV], the 
wave packet gets to the side of the billiard opposite to I 
and then it goes back to impinge again upon I. The slits 
radiates for the second time, but now the effects of the 
SC violation are evident: the coherence between %j; a and 
tp-f, is lost (see again the plot of the phase difference in 
Column III of Fig. 5). However, the probability emitted 
during this second radiation takes about 2.0 • 10 3 r to 
reach the screen. Thus, the perfect interference pattern 
in Fig. 5 corresponds to the probability piled up on the 
screen during the first radiation. This pattern then, as 
our calculations show, begins changing progressively in 
proportion to the the effects of the SC violation from 
1.3 • 10 4 r on. 

Taking up the discussion on the results in Figs. 2, 4, 
we can see that in the cases e),g), the effects of the SC 
breaking do not seem to be strong enough to lead to a 
complete decoherence of the wave packets ipa and tpb, 
whereas they do in the cases c),d),f). We illustrate this 
latter situation in Fig. 6, where we have approximated 
the intensity plots due to two incoherent waves ip a , ipb 
with the plots sum of the intensity patterns relative 
to two one-slit experiments obtained from the original 
two-slit set-up by closing alternatively the slits a and b. 
It should be recalled that we have stopped our calcula- 
tions when the probability density remained inside the 
billiards is of order 15 per cent. The further emission of 
this residual probability density could change, in some 
detail, the intensity patterns. However, we do not think 
it probable because the shapes of the intensity patterns 
shown in our figures are practically the same as those 
found when the probability leaked out from the slits was 



of order 40 per cent. 

As regards the important conclusion 2) in 1], con- 
cerning the question whether a coupling with external 
degrees of freedom is necessary or not in order that the 
quantal evolution is replaced by the classical one, our 
results show clearly that this coupling is not necessary. 
In this sense they support the theory of the "dynamical 
decoherence" proposed in On the other hand, our 
results suggest a crucial difference about the cause of 
this decoherence. Indeed, as we have seen, the wave 
packets ip a and ipj, radiated from the slits become 
incoherent not in virtue of properties depending on 
classical chaos but only in virtue of specific properties 
of the quantum-mechanical time-evolution. These 
properties are direct consequences of the violation 
of the SC (6). We can also say that our results 
show examples where the decoherence in a quantum 
system happens in virtue of properties which, besides 
being internal to the system, are fully of quantum nature. 

In conclusion, we have shown that in the class of the 
numerical experiments of interference considered in 
and in the present paper, classical chaos does not seem 
to have any role. The explanation of the results is indeed 
quite quantal: each result is determined in our experi- 
mental set-up by the degree of violation of the spatial- 
reflection symmetry w.r.t. to the y-axis of the quantum- 
mechanical time-evolution inside the billiard. This sym- 
metry violation, in turn, is determined by the analogous 
symmetry violations of the billiard shape and of the ini- 
tial Gaussian wave packet. Contrarily to the literature, 
our explanation of results is thus, in last analysis, of a 
geometric kind and it does not depend on the semiclas- 
sical nature of the system. In our numerical simulation 
we have used the semiclassical value ||fc|| = 180 just for 
intentional conformity to the work pj , but if our simula- 
tion were repeated with not semiclassical values of ||fc| (of 
course employing an experimental set-up modified suit- 
ably in order to have a good visibility of the fringes) we 
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expect that one would find out results quite analogous 
with those shown in the present paper. 
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